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We discuss a general method for constructing nonreciprocal, cavity-based photonic devices, based on match¬ 
ing a given coherent interaction with its corresponding dissipative counterpart; our method generalizes the basic 
structure used in the theory of cascaded quantum systems, and can render an extremely wide class of interactions 
directional. In contrast to standard interference-based schemes, our approach allows directional behavior over 
a wide bandwidth. We show how it can be used to devise isolators and directional, quantum-limited amplifiers. 
We discuss in detail how this general method allows the construction of a directional, noise-free phase-sensitive 
amplifier that is not limited by any fundamental gain-bandwidth constraint. Our approach is particularly well- 
suited to implementations using superconducting microwave circuits and optomechanical systems. 
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I. INTRODUCTION 

The general desire to break time-reversal symmetry and 
reciprocity in engineered photonic structures has garnered an 
immense amount of recent interest. Recall that while time- 
reversal symmetry is only a useful notion in non-dissipative 
systems, reciprocity is more general: it is defined as the in¬ 
variance of photon transmission amplitudes under exchange 
of source and detector fD On a fundamental level, the artifi¬ 
cial breaking of time-reversal symmetry allows the realization 
of truly new photonic states, such as quantum Hall states and 
more general topological states cm. On a more practical 
level, nonreciprocal devices can enable a number of signal¬ 
processing applications and greatly simplify the construction 
of photonic networks ( 71 . 

Nonreciprocal microwave-frequency devices are also cru¬ 
cial to efforts at quantum-information processing with super¬ 
conducting circuits. Here one necessarily needs to use near 
quantum-limited amplifiers to efficiently read out qubits; non¬ 
reciprocity is crucial to ensure the qubits are protected from 
unwanted noise stemming from the amplifier. The conven¬ 
tional solution is to use circulators employing magneto-optical 
effects (Faraday rotation) to break reciprocity. These devices 
have many disadvantages: they are bulky and cannot be imple¬ 
mented on-chip (hindering scaling-up to multi-qubit systems), 
and they use large magnetic fields, which can be deleterious 
to superconducting devices. Their use also typically leads to 
insertion losses. 

A number of strategies have been developed to break reci¬ 
procity without the use of magneto-optical effects in both op¬ 
tical systems and superconducting circuits. For nonrecipro¬ 
cal photon transmission, approaches based on refractive-index 
modulation ilia have been considered, as well as strategies 
using optical non-linearity in, optomechanical interaction 
ifTTl [TTll . and interfering parametric processes CMa. Re¬ 
lated strategies where the phases of external driving fields gen¬ 
erate an artificial gauge field in a lattice or cavity array have 
also been discussed (16111911 , as have alternative methods that 
do not use modulation (2Qll22ll . Nonreciprocal quantum am¬ 
plifiers have also been developed largely in the context of su¬ 
perconducting circuits (23H26ll . They typically involve engi¬ 
neering complex interferences between parametric processes. 



FIG. 1: (a) Basic recipe for generating directionality: two cavi¬ 

ties are directly coupled to one another via a coherent Hamiltonian 
77coh, and are also each coupled to the same (non-directional) dissi¬ 
pative environment, (b) The dissipative environment in (a) mediates 
a reciprocal dissipative interaction between the two cavities. This 
can be modeled using a Lindblad master equation and dissipative 
superoperator C[z] (cf. Eq. |^). By balancing the strength of coher¬ 
ent and bath-induced dissipative interactions between the cavities, 
one can break reciprocity, (c) Schematic cascaded quantum system, 
where one cavity drives another via a waveguide supporting only a 
right-propagating mode. The effective theory used to describe such 
systems corresponds to (b). 


Understanding how to achieve such interferences can be diffi¬ 
cult, though recently a graph-theory approach was formulated 
by Ranzani et al. (25]| . 

In this work, we present a simple yet general method 
for generating nonreciprocal behavior in a photonic system, 
one that can make a variety of cavity-cavity interactions 
completely directional, including amplifying interactions (see 
Fig.[^. It employs reservoir engineering (27]| . where a struc¬ 
tured dissipative environment generates useful quantum be¬ 
havior. In our approach, the dissipative reservoir (which could 
simply be a damped auxiliary cavity mode) generates an ef¬ 
fective dissipative interaction between the modes of interest. 
Nonreciprocal behavior is then obtained by balancing this in¬ 
duced dissipative interaction against the corresponding coher¬ 
ent version of the interaction. 
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As we discuss, this simple yet powerful trick allows one to 
generate both isolators (which only allow unidirectional trans¬ 
mission), and nonreciprocal, quantum-limited phase-sensitive 
amplifiers (which have zero added noise) as well as phase¬ 
preserving amplifiers (which add the quantum-limited amount 
of noise, a half quantum at the signal frequency). While our 
approach uses a kind of interference, it is markedly differ¬ 
ent from more typical interference-based approaches, in that 
it allows perfect directional behavior over a wide range of fre¬ 
quencies. The method is also simple enough that it could be 
implemented in a wide variety of architectures; in particular, 
it is extremely well suited to implementations using supercon¬ 
ducting circuits and optomechanics. 

Our approach to nonreciprocity is intimately connected to 
the theory of cascaded quantum systems ||28l [23. This is 
an effective theory developed to describe situations where a 
nonreciprocal element is used to couple two quantum systems 
(e.g., such that the output field of one cavity drives a second 
cavity, but not vice-versa, see Fig. ^c)). We show that the 
effective interactions used in this theory have exactly the form 
described above: one balances a coherent “photon tunneling” 
interaction between the two cavities against a corresponding 
dissipative version of this interaction. We also demonstrate 
that cascaded quantum systems theory is not simply an effec¬ 
tive theory for describing nonreciprocal transmission: it also 
serves as a recipe for constructing nonreciprocal devices, one 
that can be generalized to amplifying interactions. As we dis¬ 
cuss, the needed dissipative interactions can be obtained by 
simply coupling to intermediate damped cavity modes; one 
does not need to start with an explicitly nonreciprocal reser¬ 
voir (as assumed in the derivations of Refs. 1^1291 ). 

The remainder of this paper is organized as follows. In 
Sec. |I^ we introduce our basic approach of balancing coher¬ 
ent and dissipation interactions, showing how this can be used 
to generate both nonreciprocal photon transmission as well as 
amplification. In Sec.|I^ we provide further details on each of 
these schemes, discussing simple 3-mode implementations, as 
well as issues of bandwidth, impedance matching and added 
noise. We pay particular attention to our scheme for a nonre¬ 
ciprocal cavity-based phase-sensitive amplifier. In addition to 
being nonreciprocal and quantum-limited, we show that this 
system can also be constructed so that there is no fundamen¬ 
tal gain-bandwidth limitation on its performance, and so that 
it is perfectly impedance matched at both its input and out¬ 
put (i.e., there are no unwanted reflections at either port of the 
amplifier). 


II. DIRECTIONALITY FROM DISSIPATIVE 
INTERACTIONS 


Throughout this work, we consider a generic situation 
where we have a pair of cavity modes (annihilation operators 
di, (i 2 ), each coupled to input/output waveguides; our goal is 
to engineer a nonreciprocal interaction between them, thus en¬ 
abling either nonreciprocal transmission or amplification of 
signals incident on the two modes. Our approach is sketched 
in Fig.l^b): we allow both cavities to interact with one an¬ 


other in two distinct ways. The first is via a direct, coher¬ 
ent interaction described by an interaction Hamiltonian ?icoh- 
While our approach can make a general factorizable cavity- 
cavity interaction directional, we focus here on simple bilinear 
interactions. The coherent interaction will thus be described 
by a quadratic Hamiltonian, having the general form {h=l) 

^coh = Jd[d 2 + \d\d 2 + h.c., (1) 


where J and A are in general complex. We always work in a 
rotating frame where the two cavities are effectively resonant, 
and where Ticoh is time-independent. Each of the two inter¬ 
actions in this Hamiltonian could be realized in many ways; 
for example, one could start with three modes and a generic 
three-wave mixing Hamiltonian, and then displace one of the 
modes with a coherent tone. The driven modes act as pumps; 
by a suitable choice of frequencies (i.e., at the difference and 
the sum of cavities 1 and 2 resonance frequencies), one real¬ 
izes the above Hamiltonian, with the amplitudes and phases of 
the couplings J, A being controlled by the pump modes ampli¬ 
tudes. Such an approach has been exploited recently in super¬ 
conducting circuits, using the Josephson parametric converter 
(JPC) geometry ||30H3^ . as well as in quantum optomechan¬ 
ics 04l . 

The second required interaction involves controllably cou¬ 
pling both cavities to the same dissipative reservoir (Fig.l^a)). 
Eliminating this reservoir will generate an effective dissipative 
interaction between the cavities (i.e., one that cannot be de¬ 
scribed by some direct Hamiltonian coupling). The simplest 
setting is where this reservoir is effectively Markovian, and 
hence can be described using dissipators in a Lindblad mas¬ 
ter equation for the reduced density matrix p of the two cavity 
modes. As we are focusing here on a bilinear coherent interac¬ 
tion, the needed interactions between the engineered reservoir 
and the two cavities will also be linear. We are thus left with 
the general master equation 



^COh5 P 


+ r£[z]p+ 

i=i,2 


( 2 ) 


where 


2: = 


i=i,2 


( 3 ) 


and the standard dissipative superoperator C[o\ is defined as 
C[d]p =dpd^ dp - ^po'^o. (4) 

The first term in Eq. Q describes the coherent interaction 
between the two cavities, the second the interaction with the 
engineered reservoir at rate F (including the induced dissi¬ 
pative cavity-cavity interactions), and the last the coupling 
of the cavities to their input-output ports at rate nj . Note 
that an asymmetry in the couplings does not change the basic 
physics in which we are interested; thus, for simplicity, we 
take in what follows. The coefficients Uj and 

Vj characterize the individual coupling of each cavity to the 
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engineered bath. As we see in what follows, the engineered 
reservoir need not be anything too exotic: it can simply be 
another (damped) cavity mode, or a (non-directional) trans¬ 
mission line. Also note that one does not need to be in the 
strict Markovian limit, though it makes it simpler to under¬ 
stand the physics. We discuss corrections to the Markovian 
limit in Sec. Hill 

With these ingredients in place, obtaining directionality 
involves first constructing Ticoh so that it gives the desired 
behavior (amplification or transmission), and then precisely 
balancing it with the corresponding dissipative interaction 
(i.e., choice of T^Uj and Vj). To illustrate this, we can de¬ 
rive the equations of motion for the expectation values of the 
mode’s operators. Starting from the Lindblad master equation 
in Eq. ^ we obtain 


d 

dt 



d 

dt 
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with = T{\un\‘^ — \vn\‘^), {n G 1, 2) describing the local 
damping induced by the engineered reservoir, and the defini¬ 
tions /i = V 1 V 2 — U 2 u\ and u = V 1 U 2 — V 2 u\. The engineered 
reservoir mediates a non-local damping force on each mode, 
thus it couples the two modes in a similar manner as the coher¬ 
ent interaction. Crucially, due to the difference in the coupling 
coefficients we can decouple cavity 1 from cavity 2 by setting 

I r I r 

J =-111-^ X=-ijy-, (6) 

For this case we obtain a uni-directional interaction, where 
cavity 2 is driven by cavity 1 but not vice versa. Moreover, 
it is straightforward to show that this decoupling works for 
all operators: the evolution of any cavity-1 operator is inde¬ 
pendent of cavity 2, while cavity-2 operator expectations are 
infiuenced by cavity 1 (cf. Appendix [A|. 

In what follows, we show how this general recipe of bal¬ 
ancing coherent and dissipative interactions can be used to 
construct an isolator and nonreciprocal quantum-limited am¬ 
plifiers (both phase-preserving and phase-sensitive). The ba¬ 
sic recipe here will in fact allow any factorizable cavity-cavity 
interaction to become directional, including nonlinear interac¬ 
tions (see Appendix]^. It thus represents a powerful approach 
for constructing a wide variety of nonreciprocal behaviors. 


A. Unidirectional photon hopping: dissipative isolator 

We first discuss how our basic recipe can be used to obtain 
directional transmission between ports 1 and 2. We want an 
effective interaction between the two cavities which only al¬ 
lows photons to tunnel from cavity 1 to 2 (and not vice-versa). 


This is precisely the kind of behavior described by standard 
cascaded quantum systems theory |[28l|29l[35l. We show here 
how this fits into our general framework where directionality 
results from balancing coherent and dissipative interactions. 
We also show how it can be simply realized using an auxil¬ 
iary cavity or reciprocal transmission line, and thus does not 
require an explicitly directional reservoir. 

To obtain nonreciprocal tunneling between the cavities, we 
first need to identify coherent and dissipative versions of a tun¬ 
neling interaction. The coherent version is simple: choosing 
A = 0 in Eq. we obtain a standard hopping (or beam¬ 
splitter) Hamiltonian, 

^coh + h.c. = Tihop- (2) 


For the dissipative version of a hopping interaction, we 
need a zero-temperature engineered reservoir that is able to 
absorb quanta from either cavity; crucially there needs to be 
coherence between absorption of a photon from cavity 1 ver¬ 
sus cavity 2. The jump operator i in our master equation 
Eq. Q thus needs to take the form 

i ^di -h e*^d2 = 4op- (8) 


The general master equation of Eq. ^ thus reduces to 


—p = — i 

di 




op: / 


• VC\d\ H- e^^d2]p + n E 

iei,2 

( 9 ) 


The second term describes the dissipative hopping interaction: 
the engineered reservoir can absorb a photon from either cav¬ 
ity 1 or cavity 2, and there is coherence between these possi¬ 
bilities (relative phase (p). The rate for this process is T. Note 
that via a gauge transformation, the phase (p can be shifted 
into the phase of J. We thus set (p = 0 in what follows, but 
keep J complex. 

Before discussing how to engineer such a non-local dissipa- 
tor, let us discuss the consequences. Using Eq. the equa¬ 
tions of motion for mode expectation values are 
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( 10 ) 


Note that the engineered non-local dissipation in Eq. ([^ cou¬ 
ples the two cavity lowering operators in an analogous manner 
to the coherent tunneling interaction J. On a heuristic level, 
this is because the engineered reservoir gives rise to non-local 
damping: the damping force on cavity 1 depends on the ampli¬ 
tude of cavity 2 (and vice-versa). If we only have the coherent 
hopping interaction (i.e., T = 0), or only have the dissipative 
interaction (i.e., J = 0), the coupling between the cavities 
would be reciprocal. Note however that the coherent coupling 
involves J in the first line of Eq. ( p^ , and J* in the second 
line. The possibility thus emerges to have the two coupling 
terms cancel in one of the two equations. By setting, e.g.. 


7 ' r 


( 11 ) 
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Coupling \J\/T Frequency lj/hi 


FIG. 2: (a) Realization of the engineered reservoir via an auxiliary 
cavity mode that is damped at a rate n . For strong damping n ^ n 
this setup corresponds to a Markovian reservoir, (b) Implementation 
based on a transmission line, which supports propagation of photons 
in both directions, (c) Scattering matrix elements for the dissipative 
isolator setup at zero frequency, as a function of the coherent hop¬ 
ping J; the phase of J is fixed so that arg(J/r) — 'k j2. When J 
is tuned as per Eq. |TT] ), the system only allows directional transmis¬ 
sion between cavities 1 and 2. We have fixed the dissipative cou¬ 
pling strength F to be equal to the cavity damping rate and have 
taken the Markovian limit for the engineered reservoir ac). (d) 

Scattering matrix elements as a function of frequency, when the di¬ 
rectionality condition of Eq. |TT] ) is fulfilled. In the Markovian limit, 
directionality holds over all frequencies. 


we obtain a unidirectional interaction: cavity 2 is driven by 
cavity 1, but not vice-versa (see Fig.j^c)). 

With this tuning of J, our master equation Eq. takes the 
standard form used in cascaded quantum systems theory l(35]| : 



(r + «) y] L 

n=l,2 





We are most interested in the evolution of the extra-cavity 
fields, i.e., signals entering and leaving the two cavities via 
the coupling waveguides. Treating our engineered dissipative 
reservoir using a Markovian oscillator bath (which is equiva¬ 
lent to the above Lindblad description), one can use standard 
input-output theory to calculate the relation between the input 


fields incident on the two cavities, ti ngn, an d the output fields 
leaving the cavities, (in,out (see Sec. Ill A for details). Using 
the input-output boundary condition (in,out = <in,in + 


(3611371, and letting D[co’] = \ di[uj\, d 2 [oo^ , scattering be¬ 
tween the cavity in/out fields is described by a 2 x 2 scattering 
matrix s[cc;]: 


(i.e., on-resonance) scattering matrix is 


s[0] = 


T-k 

Av-j-F 

4Avr 

(/^+r)2 


0 


r-K, 

At-i-r 


(14) 


As expected, there is transmission from port 1 to port 2, but 
not vice versa. Note that s is in general not unitary, and hence 
the noise ^ must be non-vanishing in order to preserve canon¬ 
ical commutators of the output fields; we discuss this noise 
in more detail in Sec. [Ill A[ showing that it can indeed have 
the minimal amount required by quantum mechanics. We also 
show that the vanishing of Si 2 can be made to extend up to fre¬ 
quencies comparable to the relaxation rate of the engineered 
reservoir (i.e., much larger than k), see Fig.j^d). 

Equation ( p^ still does not have the ideal scattering matrix 
of an isolator la, as the incident signal on cavity 1 could be 
partially reflected. To suppress such reflections, we simply 
impedance match the system, i.e., tune T = k,. We then obtain 
the ideal isolator scattering matrix (on resonance) 


s[0] = 



(15) 


On a physical level, interference causes signals incident on 
cavity 2 to be perfectly dumped into the dissipative reser¬ 
voir. Interference also ensures that signals incident on cavity 1 
never end up in this reservoir, but instead emerge from cavity 

2 . 

We still have not specifled how one obtains the required 
non-local dissipator. The original works on cascaded quan¬ 
tum systems assumed an inherently nonreciprocal, unidirec¬ 
tional reservoir (i.e., a chiral transmission line), and then de¬ 
rived the effective master equation of Eq. ( p^ . However, the 
above dynamics can be obtained without needing an explic¬ 
itly directional reservoir. One simple choice would be a one¬ 
dimensional transmission line, cf. Eig. [^b), supporting both 
right-moving and left-moving modes, which couples to cavity 
j at position Xj 

AisB = -\j^^ X] + ^Lixj)] + /i.c.) , 

* i=i,2 

(16) 

where cl^r{x) denote the left and right moving flelds in the 
waveguide, and vq is the waveguide velocity. Eor a suitable 
choice of \xI — X 2 \, one again obtains the master equation of 
Eq. ^ with jump operator z as per Eq. Eurther details are 
provided in the Appendix]^ 

Another simple implementation involves taking a damped 
auxiliary mode as the engineered reservoir (annihilation oper¬ 
ator c), see Eig. [^a), which interacts with the two principle 
modes via a Hamiltonian: 

'Hsb = J'c^ (di + d 2 ^ + h.c.. (17) 


DoutM = sM DinM + (13) 

Here, ^[cc] describes (operator-valued) noise incident on the 
cavities from the engineered reservoir, and the zero frequency 


Such a quadratic interaction can be realized in a tunable fash¬ 
ion by starting with a three-wave mixing Hamiltonian and 
pumping one of the modes at an appropriate frequency; this is 
the same strategy used to implement the coherent direct inter¬ 
action in Eq. 0 (see discussion following that equation). As 
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we show in Sec. Ill A if the damping of the auxiliary mode k' 
is sufficiently large, it can be adiabatically eliminated, yield¬ 
ing the scattering matrix given above. 

For this particular realization, our isolator reduces to a 
three mode system with an asymmetric choice of damping 
rates. Furthermore, the required phase of J in the direction¬ 
ality matching condition of Eq. ([nj corresponds to having 
the three mode system pierced by an effective magnetic flux 
of a quarter fiux quantum. We stress that in many physical 
implementations, the couplings J and J' are tunable simply 
by controlling the amplitude and phases of the relevant pump 
modes. Thus, the directional interaction we finally obtain is 
not the result of having used an explicitly directional reser¬ 
voir, but rather it results from the control of relative phases in 
a driven system. Note that this three-mode realization of our 
dissipative isolator was also discussed by Ranzani et al. ||25l. 
It is also interesting that this three-mode implementation di¬ 
rectly yields the scattering matrix of an ideal circulator (see 
Sec. Ill A ); it is closely analogous to previous proposals for 
non-magnetic circulators ll^l^[38llT9ll . 






Frequency uj/[n{l — C)] 


Frequency oj/[k{1 — C)] 


Cn 

IjT 


FIG. 3: s-matrix elements of the directional, phase preserving am¬ 
plifier, as a function of scaled frequency; the cooperativity C = 
T/k = 0.95, where F is the dissipative interaction strength, and 
K, is the damping rate of cavities 1 and 2. (a) Auxiliary mode damp¬ 
ing K = K, indicating a strong deviation from the Markovian limit; 
while perfect isolation exists at a; = 0, it is rapidly lost for non-zero 
frequencies, (b) Auxiliary mode damping k' = lOOn, closer to the 
Markovian limit. The directionality is much better at finite frequen¬ 
cies, while the gain is unchanged. 


B. Directional phase-preserving quantum amplifier 

We next use our general recipe to construct a nonrecipro¬ 
cal, phase-preserving amplifier, a topic that is of considerable 
interest to the superconducting qubit community ||23H25]| . We 
again consider a two-mode system as sketched in Fig. [^b). 
Our goal is a dynamics that leads to signals incident on cavity 
1 emerging amplified from cavity 2, while at the same time, 
signals (and noise) incident on cavity 2 are prevented from 
emerging from cavity 1. 

Our basic recipe is the same as the previous subsection: en¬ 
gineer both coherent and dissipative versions of the desired in¬ 
teraction, and then balance them to obtain directionality. The 
coherent interaction needed corresponds to a non-degenerate 
parametric amplifier (NDPA), as obtained by setting J = 0 in 
Eq.Q: 


Hcoh ^Xd\dl + X*dj2 = Hpa. (18) 

This textbook interaction results in the amplification of an in¬ 
put signal (or noise) incident on either cavity, in both trans¬ 
mission and refiection (see, e.g., Refs. [J7ll40l) . 

We next need to add the dissipative version of this NDPA 
interaction, as mediated by an appropriately chosen dissipa¬ 
tive reservoir. This kind of dissipative amplification was re¬ 
cently introduced in our previous work. Ref. BTl . The dissi¬ 
pative reservoir now needs to be able to absorb photons from 
one cavity and to emit photons to the other, with coherence 
between these possibilities. The jump operator z associated 
with the reservoir (cf. Eq. Q) thus needs to take the general 
form 


z -^V2 (^cosOdi + sin = zpa, (19) 

where the angle 0 parametrizes the asymmetry between the 
two kinds of processes. The relative phase cp can again be 


gauged away into the phase of A; we thus set it to zero in what 
follows. 

With this choice of coherent Hamiltonian and dissipator, the 
two-cavity system is again described by the master equation 
Eq. 0., with F parametrizing the strength of the coupling to 
the engineered reservoir, and hence of the dissipative ampli¬ 
fier interaction. To see clearly that the dissipation here leads 
to amplification, we consider the equations of motion for the 
means of lowering operators. One obtains 


dt 

d 

dt 



+ 2F cos^ ^ / i \ 

F 

— sm20 + iX 

( 4 ), 

2 \ V 

K-2rsin2 0 

^ sin 2^ + iA* 

] 

2 \ V + 



(20) 


The crucial terms behind the amplification are the last term in 
each line, which cause di to be driven by d^ and vice-versa. 
Again, both the coherent interaction and the dissipative inter¬ 
action give rise to such terms; each interaction thus facilitates 
amplification that can be quantum limited, but that is not di¬ 
rectional ED. To obtain a unidirectional interaction, we again 
simply tune the amplitude and phase of the coherent interac¬ 
tion with respect to the dissipative interaction, so as to cancel 
the coupling term in the first equation, i.e., 

X = i— sm20. (21) 

2 

To see that this choice gives the desired behavior of the 
output fields, we model the dissipative bath as a Markov 
reservoir, and calculate the scattering matrix for the sys¬ 
tem using input-output theory. Letting Din/out [^] = 

>4, in/out ’ ^^6 input-output relations take 

the form of Eq. (0. I* again describes noise incident from 
the engineered reservoir, while the 2 x 2 scattering matrix s 
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takes the following explicitly nonreciprocal form at zero fre¬ 
quency 


s[0] = 


2C cos^ e-1 

2C cos^ 9-\-l 

_ AC sin 29 _ 

[2C cos2 0+l][2Csin2 9-1] 


" V 

2Csin^ 9-\-l / 

2Csin2 0-l / 


Here, the cooperativity is given as C = ^. 
If we further tune 0 so that 


( 22 ) 



cos^e = 1/{2C), (23) 

(possible as long as C > 1/2), we cancel all reflections of 
input signals incident on cavity 1. With this tuning, the scat¬ 
tering matrix becomes 


FIG. 4: Schematic illustrating the directional phase-sensitive ampli¬ 
fier. The coherent QND Hamiltonian of Eq. causes the cavity-1 
P quadrature to drive the cavity-2 P quadrature, and the cavity-2 X 
quadrature to drive the cavity-1 X quadrature (blue arrows); there is 
gain associated with each of these drivings, cf. Eqs. ( [Z7] ). The engi¬ 
neered reservoir (jump operator described by Eq. ( [26] )) also mediates 
the same drivings (green and magenta arrows), again with gain. By 
balancing these interactions, one can cancel the X 2 Xi driving, 
resulting in directional amplification. 


with Q = • As claimed, we have a scattering matrix de¬ 

scribing nonreciprocal, phase-preserving amplification, with 
a gain that diverges as C approaches 1. Signals incident on 
cavity 1 are never reflected, and emerge from cavity 2 with an 
amplitude gain S 21 = VG, whereas signals incident on cav¬ 
ity 2 do not emerge at the output from cavity 1. The system 
exhibits a standard parametric instability when C > 1 (analo¬ 
gous to the instability in a standard, coherent NDPA). 

The frequency dependence of the scattering coefficients is 
discussed in Sec. |IIIB| Strikingly, the directionality property 
si 2 [<^] = 0 holds for all frequencies for which the Marko¬ 
vian bath approximation is valid. The system is limited by a 
standard gain-bandwidth constraint (in contrast to the purely 
dissipative amplification process, which is has no such con¬ 
straint ED). We also discuss the added noise of the ampli¬ 
fier in Sec. |IIIB| showing that it is quantum limited in the 
large gain limit as long as there is no thermal noise incident 
on cavity 2; surprisingly, the engineered reservoir need not be 
at zero-temperature. 

While there are many ways to realize the engineered reser¬ 
voir used in this scheme, the simplest choice is a damped third 
auxiliary mode, see Sec. |IIIBl With this particular choice, our 
scheme reduces to the 3-cavity amplifier discussed by Ranzani 
and Aumentado in Ref. |25| Our analysis thus generalizes this 
scheme, and provides insight into the underlying mechanism. 
It also shows the crucial importance of having the auxiliary 
mode damping hz' be much larger than that of the principle 
modes; in this Markovian limit, one has directionality over 
the full amplification bandwidth (see Fig.[^. 


C. Directional phase-sensitive amplifier 

As a third application of our recipe for nonreciprocity, we 
construct a phase-sensitive amplifier. Phase-sensitive ampli¬ 
fiers only measure a single quadrature of an incident signal; 
as a result, quantum mechanics allows them to amplify with¬ 
out adding any added noise Eiiiia. Our general approach 
allows one to construct a nonreciprocal and noiseless version 
of such an amplifier, again using the two-cavity-plus-reservoir 


setup in Fig. [Jb). The resulting amplifier has another strik¬ 
ing advantage over a standard paramp: it does not suffer from 
any fundamental gain-bandwidth limitation, a point we dis¬ 


cuss more fully in Sec. Ill C 


As before, the first step is to construct a coherent interac¬ 
tion that gives the desired amplification. The standard choice 
would be a degenerate parametric amplifier (DPA) Hamil¬ 
tonian involving just a single mode, of the form Hint = 
Xdd + A*d!^d!^. In contrast, to be able to implement our recipe 
for directionality, we want an interaction that couples two 
modes. 

Surprisingly, there is a simple coherent two-mode interac¬ 
tion which does the job and which yields ideal amplification 
properties (zero added noise, no gain-bandwidth limitation). 
One needs to use the kind of quantum non-demolition (QND) 
interaction discussed extensively in the context of Gaussian 
cluster-state generation |[43] - l45l . Suppose we want to am¬ 
plify the P quadrature of cavity 1, i.e., the operator Pi = 


—— d\)/y^. We then use an interaction Hamiltonian 
that commutes with this operator, but that takes information 
in the Pi quadrature and dumps it into a cavity 2 quadrature 
(i.e., P 2 ). The required coherent Hamiltonian is obtained by 
setting J = A = iAQND/2 (Aqnd ^ in Eq. i.e.. 


Hcoh ^Aqnd P 1 X 2 = Hqnd, (25) 


with X2 = ^(^2 + ^2 j /v^- It is straightforward to see from 

the Heisenberg equations of motion that Hqnd causes P 2 to 
be driven by Pi , and hence P 2 will contain information on Pi 
(see Fig.|^. The same holds for the extra-cavity fields: the P 
quadrature of a signal incident on cavity 1 will emerge in the 
P quadrature of a signal leaving cavity 2. 

Note that Pi and X 2 are QND variables: they commute 
with the Hamiltonian in Eq. ( [25] ), and are undisturbed by the 
amplification process. It follows that there is no possibility 
of feedback in this system, and hence the system is stable ir¬ 
respective of the value of Aqnd- By increasing Aqnd, one 
can thus achieve increasing amounts of phase-sensitive gain. 
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Furthermore, as the amplification mechanism here does not 
involve coming close to an instability, the amplification band¬ 
width is always ^ k, irrespective of the gain. 

Following our recipe for directionality, we next need to con¬ 
struct the dissipative counterpart to the coherent interaction in 
Eq. ^25\ . We need the jump operator z characterizing the en¬ 
gineered reservoir to also preserve the QND structure of the 
coherent Hamiltonian. Taking m = U 2 = V 2 = \/2 and 
vi = — \/2 in Eq. ^ yields 

Z -^X2 + iPl = ^QND- (26) 


This dissipative interaction is the counterpart of the coher¬ 
ent interaction in Eq. ( [^ : with this choice of i, the dissi¬ 
pative terms in Eq. Q alone lead to amplification of the P 
quadrature of signals incident on cavity 1. The heuristic in¬ 
terpretation of this dissipative amplification is similar to that 
presented in Ref.|4T]for the phase-preserving case: the engi¬ 
neered reservoir “measures” the QND quadrature Pi, and then 
dumps this information into the non-QND quadrature P 2 (see 
Fig.|^ as well as Sec. |III Q . 


With these choices for Pcoh and i in Eq. we have both 
coherent and dissipative phase-sensitive amplifying interac¬ 
tions between the cavities. Using this master equation, the 
equations of motions for the quadrature means have the ex- 


pected form: 



d 

/ - \ 

n 


dt 

(^0 


(ft), 

d 


n 


dt ^ 

(x,) 


(x,). 

d 


n 


dt ^ 

(x.) 


(x.) + 

d 


n 


dt 

( ft ) 


(ft,)- 


(27) 


Pi and X 2 are QND variables and thus undisturbed by either 
interaction. In contrast, both interactions cause P 2 to become 
an amplified copy of Pi. 

We can now finally apply the last step of our general recipe: 
balance the dissipative and coherent interactions to break reci¬ 
procity. This simply involves setting 


r = A 


QND, 


(28) 


which ensures that cavity 1 is insensitive to the state of cavity 

2 . 

Finally, we are as usual interested in the behavior of the out¬ 
put fields from the cavity. Treating the engineered reservoir as 
a Markovian bath and using input-output theory, we can again 
calculate the scattering matrix of the system. Writing this ma¬ 
trix in a quadrature representation, we find that on-resonance 

(i.e., at zero-frequency) Zout = s Zin + ^ with 



FIG. 5: Gain and reverse gain of the directional phase- 
sensitive amplifier, plotted as a function of signal frequency u and 
coherent coupling strength Aqnd (cf. Eq. ([^), assuming that the 
dissipative coupling E always satisfies the matching condition F = 
Aqnd (cf. Eq. ([^). Q(i)[uj] describes the amplification in transmis¬ 
sion of signals incident on the cavity- 1 P quadrature, while G(p[uj] 
describes the amplification in transmission of signals incident on the 
cavity-2 X quadrature. We have taken the engineered reservoir to 
be an auxiliary cavity mode with damping rate n /n = 100 (see 
Sec. |III C| ). In this limit, deviations from the Markovian-reservoir ap¬ 
proximation are small. 


Here, the zero-frequency amplitude gain is given by ^/G^ = 

8Aqnd 

K 

The input-output relations in Eqs. ( [29| ) describe an ideal di¬ 
rectional degenerate amplifier: the P quadrature of signals 
incident on cavity 1 emerge with gain from cavity 2, whereas 
signals or noise incident on cavity 2 never emerge from cav¬ 
ity 1. Note further that there is no unwanted amplification in 
refiection of incident signals and noise. The amplifier also 
has several other remarkable properties: it is quantum lim¬ 
ited (i.e., no added noise in the large gain limit), and does not 
suffer from any fundamental gain-bandwidth limitation. The 
directionality is also maintained over a large range of frequen¬ 
cies (see Fig.[^. These properties (along with the possibility 
of eliminating unwanted reflections) are discussed in more de¬ 
tail in Sec. lIIICI 


III. NOISE, BANDWIDTH AND THREE-CAVITY 
IMPLEMENTATION 

A. Dissipative isolator: additional details 

7. Auxiliary-cavity implementation of the engineered reservoir 


s[0] 



-1 

0 

0 

0 

\ 

/w\ 


0 

-1 

0 

0 



A 


0 

0 

-1 

0 


, ^ — 



0 


0 

-1 



V A / 


(29) 


To demystify the engineered reservoirs used in our 
schemes, we provide more details here on the simplest pos¬ 
sible realization: a damped auxiliary cavity mode. In the limit 
where the damping rate n' of this auxiliary mode is large, this 
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model describes a general Markovian reservoir. We stress that 
this setup is just one of many ways to implement the necessary 
dissipative dynamics. In Appendix [B| we explicitly show how 
coupling two cavities to a (non-directional) one-dimensional 
transmission line or waveguide also generates the needed dis¬ 
sipative dynamics. 

Consider the dissipative isolator described by the master 
equation in Eq. (|^, and take the engineered reservoir to be an 
auxiliary mode with lowering operator c which is damped at 
rate n' by a coupling to a Markovian reservoir. 1/ n' will act as 
the correlation time of our engineered reservoir. As discussed 


in Sec. IIA we need this auxiliary mode (i.e., the engineered 
reservoir) to interact with the principle modes via the interac¬ 
tion Hamiltonian in Eq. (17). The simplest limit is where k' is 
much larger than all other frequency scales; in this limit, the 
c mode will itself act as a Markovian reservoir for the system 
modes ^ 1 ,^ 2 - One could then recover the master equation of 
Eq. ([^ using standard adiabatic elimination techniques Ea. 

Alternatively, one can eliminate the auxiliary mode within 
a Heisenberg-Langevin formalism, using the coherent Hamil¬ 
tonian i-L = Hhop + ^SB • Solving the equation of motion for 
c in the large-damping (adiabatic) limit yields 


= Cin Ml +^2 , 


(30) 


where all operators are evaluated at the same time, and cin 
describes thermal and vacuum fluctuations stemming from the 
mode’s internal dissipation. Substituting this equation into the 
equations of motion for the principle cavity operators dj then 
yields: 


d c- /i: “h r ^ !— '1 . 

— ^ d\ V in “1“ '2'V ECin 

d ^ “t“ r ^ !— '1 . /t^/^ 

~U2 = ~ U2 v ^U2 in iv ECin 

dt 2 


-|- iJ 


- + .r 


di, 

(31) 


where we take J' G M without loss of generality, and deflne 
r = Taking average values, we recover the master- 

equation result of Eq. ( p^ . 

Using the Heisenberg-Langevin approach, we can now cal¬ 
culate ihc full scattering matrix for the system which includes 
the scattering of noise incident from the engineered reservoir. 

Letting Y[ijj] = ^di[cc;], d 2 [^], ^ the full scattering re¬ 

lations take the form Yout[^] = ^[cc] Yin[co’]. Consider first 
the Markovian limit, where k' ^ uj^k^T. Assuming that the 
system has been tuned to satisfy both the directionality condi¬ 


tion J = (cf. Eq. \ll\) and the impedance matching con¬ 


dition /^ = r, the full scattering matrix in this limit is 




s w = 


— 2 — 

1 






\ 




l-i^ 


1-i^ 

uw 

(^) 


+0 


(32) 


The upper left 2x2 matrix is the scattering matrix s for the 


reduced, two-mode system, cf. Eq.(p^. The elements 513 and 
523 describe the scattering of noise from the engineered reser¬ 
voir to the main cavity modes. This then explicitly yields the 


noise operator in Eq, 


■ 0 


as ^ = [si 3 , § 23 ]’^Cin. We see that 


directionality holds for all frequencies in this Markovian limit, 
i.e., 5 i 2 [co’] = 0. In contrast, the impedance matching (which 
ensures no reflections at the input of cavity 1 ) only holds for 

uj n. 

Einally, note that at zero frequency, the full scattering ma¬ 
trix becomes: 


s[0] = 


0 

0 

A 

1 

0 

0 

, 0 

i 

0 . 


(33) 


In this ideal case, signals incident on cavity 2 are perfectly 
transmitted to the engineered reservoir, while the input held 
on the reservoir (i.e., the c mode) is perfectly transmitted to 
mode 1. If the engineered reservoir is at zero temperature, 
we see that the output from cavity 1 is simply vacuum noise. 
Amusingly, the above unitary scattering matrix is that of a per¬ 
fect circulator: the effective magnetic held associated with the 
phase of J breaks the degeneracy of right and left circulating 
eigenmodes of the coherent three-mode hopping Hamiltonian. 
In the case of symmetric decay rates, i.e., this kind of 

circulator has been discussed in the context of superconduct¬ 
ing circuit setups Eoiiaa and just recently experimentally 
demonstrated by Sliwa and co-workers 1^ . An analogous 
circulator for phonons has been discussed in the context of 
optomechanics 


2. Non-Markovian corrections 

We can also consider deviations from the Markovian limit, 
where the internal damping rate of the engineered reservoir 
tv' is not arbitrarily large. The scattering matrix follows sim¬ 
ply from solving the full (linear) Langevin equations without 
any adiabatic assumption. We quote only the results for the 
forward and reverse transmission probabilities, again assum¬ 
ing that the directionality and impedance matching conditions 
have been met. We And 


=T 




4(i+^)-ls:^, + (i+S)' 


I«5i2[c^J| =t 


2x2- 


/i _i_ 4^\ _ j_ 1 _L 

) 


(34) 


One clearly sees that the directionality only holds for frequen¬ 
cies that are small compared to the inverse correlation time 
1/tz' of the reservoir: for small uj, |5i2[co’]p (X uj‘^ltv'‘^. Eor 
non-zero uj/k,' , the engineered reservoir gives rise to both dis¬ 
sipative and coherent interactions. The extra induced coher- 
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ent interaction ruins the directionality matching condition of 
Eq. ([nj, leading to a lack of perfect isolation. 


B. Directional phase-preserving amplifier: additional details 

1. Bandwidth and non-Markovian effects 


We return now to the setup presented for directional am¬ 
plification in Sec. |IIB[ As in the previous section, we will 
investigate the frequency-dependent behavior of the system 
using an auxiliary damped cavity mode c to represent the en¬ 
gineered reservoir. With this choice, the system is analogous 
to that studied in Ref. [25], which was recently implemented in 
a superconducting circuit experiment 1^ . We emphasize the 
importance of having a large damping rate k,' of the auxiliary 
mode, thus complementing the discussion in Ref.[25l 

For phase-preserving amplification, the coherent interac¬ 
tion between the principle modes di , ^2 has the NDPA form of 
Eq. ( p^ . To obtain the correct dissipative interaction, the cou¬ 
pling HsB between the principle cavity modes and the auxil¬ 
iary mode should have the form 

TisB = V2X'c^ (^cosOdi ffsmOdl^ -\-h.c.. (35) 


Taking the large tz' limit and using standard adiabatic elimi¬ 
nation techniques, one recovers the master equation described 
by Eqs. ([^ and (19), with F = 

One can again solve the full Heisenberg-Langevin equa¬ 
tions to obtain the full 3x3 scattering matrix for the system at 
all frequencies. If we tune the couplings to satisfy the direc¬ 
tionality condition of Eq. ( [21] ) and the impedance matching 
condition of Eq. the “forward photon number gain” of 
the amplifier takes the form 


Q[w\ = |S2iM|^ 



(2C - 1) 
\c-if + 



+ 0 



(36) 


The corresponding reverse photon number gain (which we 
ideally want to vanish) is given by 


with the price that this interaction naturally leads to a conven¬ 
tional gain-bandwidth limit. 

Turning to the non-Markovian effects, we see from Eq. 
that for finite u;/k\ the reverse gain is non-zero, implying that 
directionality is lost; this is also depicted in Eig.[^ The loss of 
directionality here is analogous to what happens in the direc¬ 
tional isolator, and occurs for the same basic physical reason: 
for finite cc/^c', the engineered reservoir also induces a coher¬ 
ent interaction between the two modes, and hence the perfect 
matching of coherent and dissipative interactions needed for 
directionality is lost. 


2. Added noise and quantum-limited behavior 


In addition to directionality, for many applications it is cru¬ 
cial that our amplifier reaches the fundamental quantum limit 
on its added noise. This limit corresponds to adding noise 
equivalent to half a quanta at the input, fiadd > 1/2 14^ . 
The added noise follows directly from the full scattering ma¬ 
trix 5, and will have contributions both from noise incident on 
cavity 2 that is reflected, and noise emerging from the engi¬ 
neered reservoir (i.e., the auxiliary c mode). Assuming that 
the impedance matching and directionality conditions have 
been fulfilled, and letting n and represent the thermal oc¬ 
cupancies (respectively) of these two noise sources, we find: 


'^add [0] — 




(38) 


Thus, in the large gain limit, the added noise is quantum lim¬ 
ited as long as there is no thermal noise incident upon cav¬ 
ity 2 (i.e., = 0) (251. Remarkably, thermal noise in the 

engineered reservoir does not prevent one from reaching the 
quantum limit; similar behavior is found in a purely dissipa¬ 
tive (non-directional) phase-preserving amplifier ED. 

While not relevant to the quantum limit, from a practical 
standpoint one also wants the noise leaving cavity 1 to be 
small (so as not to damage the signal source). Using our scat¬ 
tering matrix, it is straightforward to calculate the noise of the 
output field from cavity one. Characterizing this noise by an 
effective thermal occupancy out’ we find at zero frequency: 


G[0J] = |S12HP 




(37) 


Consider first the limit where the engineered reservoir is ef¬ 
fectively Markovian, uj/k,' 0. The reverse gain always van¬ 
ishes, while the zero frequency forward gain ^[0] is controlled 
by C, and diverges as C 1; the system is unstable for larger 
C. In the large gain limit, Q[uj] is a Lorentzian as a function of 
frequency, with a bandwidth Acc = 2k,{1 — C) that decreases 
as one increases the gain.. The amplifier has a finite gain- 
bandwidth limitation just like a standard cavity-based NDPA 
(i.e., the product is fixed) f25\ . Note that the dissi¬ 

pative parametric interaction on its own suffers from no such 
limitation ED, but is of course not directional. Directionality 
is thus obtained by introducing a coherent NDPA interaction. 


• (39) 

Thus, while thermal noise in the engineered reservoir does 
spoil quantum limited performance, this noise does show up 
in the output of cavity 1. 


C. Directional phase sensitive amplifier: additional details 

1. Full scattering matrix and impedance matching 

We now turn attention to our scheme of Sec. Ill Cl for direc¬ 
tional and noiseless single-quadrature amplification. As dis¬ 
cussed in that section, we need to combine the coherent QND 
interaction of Eq. ( [25] ) (QND variables X 2 and Pi) with the 















10 


corresponding dissipative interaction; this dissipative interac¬ 
tion requires the jump operator z = X 2 -h iPi, as given in 

Eq. 

To generate the required dissipation, we again take the en¬ 
gineered reservoir to be a damped auxiliary mode c (damp¬ 
ing rate k,'). Writing this operator in terms of quadratures as 
c = (t/ + iV)ly/2, the required system-bath interaction has 
the form 


= A 


PiV^X2U 


(40) 


The general structure of Eq. ( [4T] ) implies that once the aux¬ 
iliary mode is eliminated, the cavity Xj modes will drive one 
another in a non-directional way; the same goes for the cavity 
Pj quadratures. This is depicted schematically in Fig. |^a). 
To obtain directionality, we need to cancel the ability of 
the cavity-2 quadratures to drive the corresponding cavity-1 
quadratures. We do this in the usual manner: we balance 
the dissipative quadrature-quadrature interactions generated 
by Eq. ( [4T] ) by the coherent versions of these interactions. This 
will require a coherent Hamiltonian of the form: 


This interaction preserves the QND structure in the coherent 
interaction, as it also commutes with X 2 and Pi- One can 
again confirm that in the Markovian limit of a large k,', one 
recovers the master equation description, with the dissipative 
rate T in Eq. (|^ being given by T = 

The dissipative interaction on its own generates phase- 
sensitive amplification that can be quantum limited. Heuristi- 
cally, this can be understood as arising from a kind of trans¬ 
duction mediated by the reservoir. From Eq. ( [40| information 
in the Pi quadrature of cavity 1 drives the auxiliary mode U 
quadrature. The U quadrature in turn drives the cavity-2 P 2 
quadrature, effectively letting Pi drive P 2 . As Pi and X 2 
are QND variables, increasing A simply increases the gain as¬ 
sociated with this process, with no possibility of instability. 
An analogous argument of course shows that one obtains re¬ 
verse gain: signals incident on X 2 will emerge amplified in 
Xi. Thus, the dissipative amplification here is not directional; 
directionality is only obtained when this process is matched 
against its coherent counterpart, as described in Sec. |IIC[ 

The only small non-ideality left in the directional phase- 
sensitive amplifier of Sec. |IIC| is the presence of refiections 
at the input (cf. Eq. ([29|)). Even though there is no gain asso¬ 
ciated with these, one would ideally want them to be exactly 
zero to protect the signal source. We now show that this can 
be easily accomplished by modifying both the coherent and 
dissipative interactions used in the scheme, so as to slightly 
deviate the QND structure discussed above. This modifica¬ 
tion also allows one to cancel refiections of signals and noise 
at the amplifier output. 

To impedance match, we first modify the system-bath 
Hamiltonian in Eq. ( [40| to take the more general form 

PsB = V2Ku U ^sin^Ai -f cos ^^ 2 ^ 

+ V2Av V (^cosOPi + sin6>P2^ . (41) 


For ^ = 0 and Au = Ay = A/\/2 we recover Eq. (40). By 
allowing Ajj ^ Ay , we modify the relative strength of the 
two QND interactions. By letting d deviate slightly from zero, 
we break the QND structure of Eq. ( [40| ); this will allow us to 
cancel the unwanted refiections from both cavities. In what 
follows, it will be useful to parametrize the system-bath cou¬ 
plings in terms of a cooperativity C and asymmetry parameter 
a\ 


C = 4At7Ay/(/^/i:'), (42) 

a = {AvlAup, (43) 


Pcoh =\lPlX 2 + A2P2-A1. 


(44) 


The second term here is new compared to Eq. ( [25] ), and breaks 
its QND-structure. As usual, we balance the above coher¬ 
ent interactions against their dissipative counterparts (as gen¬ 
erated by Eq. ( |4T] )) so that the cavity-1 quadratures are not 
driven by the cavity-2 quadratures. Working through the equa¬ 
tions of motion, and focusing on the Markovian limit, we ob¬ 
tain the directionality conditions 


Ai = nC cos^ 6>, A 2 = —tiC sin^ 0. 


(45) 


Using a standard Heisenberg-Langevin analysis, we find the 
full scattering matrix of the system; tuning the coherent inter¬ 
actions as per Eqs. (^), the scattering will indeed be direc¬ 
tional. Insisting further that there are no refiections of signals 
and noise incident on either cavity (i.e., impedance matching) 
leads to an additional condition on the angle 0: 


sm20 = 1/C. 


(46) 


Note that for a large cooperativity C, the angle 0 is very close 
to zero, implying that one is very close to our original scheme 
where Pi and X 2 are QND variables. 

Even after satisfying the above conditions, the cooperativ¬ 
ity C and the asymmetry parameter a remain unspecified; they 
control the final form of the impedance-matched, directional 
scattering matrix. Using the above conditions, the full scatter¬ 
ing matrix of the system (describing both the principle modes 
and the auxiliary mode c) takes a simple form. Introducing 
the vector 


W = {Xi Pi X 2 P 2 u vf, (47) 

the scattering relations at each frequency then take the form 
Wout = s Win. At zero frequency, the scattering matrix is 
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( 48 ) 
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FIG. 6: Properties of the directional phase-sensitive amplifier, (a) Sketch of couplings and drivings used to impedance match the amplifier, 
(b) Reverse gain Q(f)[u\ for uj set to half of the amplification bandwidth Acj, for various choices of the auxiliary-mode damping rate On 
resonance we always have perfect directionality: ^(^[0] = 0. (c) Amplification bandwidth Acj as a function of zero-frequency forward gain 
Q(j) for various values of n . The amplifier does not suffer from a standard gain-bandwidth constraint. 


Q(j) describes the zero-frequency phase-sensitive photon num¬ 
ber gain of our amplifier, and it is given by 


= C 




(49) 


1 emerging in the P quadrature of the output from cavity 2. 
Assuming that we chose parameters to impedance match (as 
described), and that we further tune the asymmetry parame¬ 
ter a to minimize noise as per Eq. ( [5Q| , the forward photon 
number gain is given by: 


The upper 4x4 block describes an ideal, directional phase¬ 
preserving amplifier. As for the full 6x6 scattering matrix, 
it describes a kind of “squeezing circulator”, where the input 
on port j emerges from port j + 1 after having undergone 
a squeezing transformation. Crucially, the squeezing parame¬ 
ters or gains for each of these transformations are not all equal. 

While a squeezing circulator may have interesting appli¬ 
cations, if the goal is amplification, it represents a potential 
hazard. As indicated by Eq. ( [48] ), incident noise from the aux¬ 
iliary mode will emerge from the output of cavity 1, having 
undergone a squeezing transformation with gain To 

protect the signal source at the amplifier input, we do not want 
to amplify any fluctuations emerging from the auxiliary mode. 
Hence, the ideal choice is to null this effective gain by tuning 
the asymmetry a to satisfy 


= l/G^, (50) 

Tuning the asymmetry of the couplings to the auxiliary cav¬ 
ity in this manner ensures that one can have large, directional 
gain, without unduly large amounts of noise emerging from 
the amplifier input port. 

The presented phase-sensitive amplifier has several highly 
desirable properties: it is quantum limited, directional and 
has no gain-bandwidth limitation. However, an experimental 
implementation in a superconducting circuit setting will also 
face some technical challenges. Most notably, a straightfor¬ 
ward implementation requires 6 pump tones to be applied with 
excellent control over their amplitudes and relative phases. 
While demanding, experiments with analogous levels of com¬ 
plexity and multiple pumps have recently been performed in 
circuit QED architectures, see e.g., Refs.[39j[46j 


= 1^42 1 ^ = 


G4> (l + 




4a;4 


(51) 


As already discussed, the zero-frequency gain can be made 
arbitrarily large by simply increasing the various couplings 
(i.e., C); the linear system never exhibits any instability. In 
the Markovian limit k.' ^ uj, the frequency dependence of the 
gain is extremely simple: it is simply a Lorentzian squared, 
with a bandwidth Acc ^ n which is independent of the zero- 
frequency gain. Thus, this amplifier is not constrained by any 
fundamental gain-bandwidth limitation. 

Including non-Markovian effects (i.e., finite ujjn'), the fre¬ 
quency dependence is slightly more complex, but the ultimate 
bandwidth is still set by n, irrespective of the size of the zero- 
frequency gain. 

While deviations from the Markovian limit do not degrade 
amplification, they impact the directionality of the amplifier. 
In the ideal Markovian limit, signals incident on cavity 2 in 
either quadrature never emerge from cavity 1. Eor finite 
this is no longer true: now, the reverse-gain scattering matrix 
element sis becomes non-zero, implying that incident signals 
on X 2 can emerge from Xi . We find: 


Q^[w\ = Isisl^ 



( 1 +^) +^(1 + 


4^\ 

) 


4uj^ 

K^K,' 


. (52) 


As expected, the optimal situation is clearly in the Markovian 
limit where n' n. In this limit, one has purely directional 
amplification over the full bandwidth n of the principle cavity 
modes. 


2. Frequency dependence 

The full scattering matrix can also be easily calculated at 
non-zero frequencies. The relevant forward gain of the am¬ 
plifier describes signals incident in the P quadrature of cavity 


3. Added noise for finite frequency 

An ideal phase-preserving amplifier can amplify a single 
quadrature without any added noise 1421 . Erom the scattering 
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relations of Eq. ( [48] ) it immediately follows that our scheme 
reaches this quantum limit on resonance. For completeness, 
we also present the added noise at finite frequency, again fo¬ 
cusing on the impedance-matched version of the amplifier. 
We calculate the added noise of our amplifier in the standard 
manner Oil, by calculating the noise in the cavity-2 output P 
quadrature (symmetrized spectral density Sp^[uj]), and then 
referring this back to the input. Expressing this added noise 
as an effective number of thermal quanta, and focusing on the 
Markovian limit k' ^ we find 


^P2,add [^] 







M±i)\ 

) 


o 

(53) 


Note that we have left the asymmetry parameter a 
(cf. Eq. ( [43] )) unspecified here. The added noise always van¬ 
ishes at cd = 0, irrespective of the gain. Furthermore, for a 
fixed value of a the added noise vanishes at all frequencies in 
the large gain limit, implying that one is quantum limited at all 
frequencies. If however one tunes a = l/^^ to minimize the 
noise hitting the input port, then the added noise is non-zero 
at finite frequencies even in the large gain limit. 


IV. CONCLUSION 


We have presented an extremely general yet simple method 
for achieving directional behavior in coupled photonic sys¬ 
tems, based on matching a given (reciprocal) coherent inter¬ 
action with the corresponding dissipative version of the inter¬ 
action. We demonstrated how this principle could be used to 
construct both isolators and directional, quantum-limited am¬ 
plifiers. In particular, our approach allows the construction of 
a directional phase-sensitive amplifier that is not limited by a 
standard gain-bandwidth constraint. The recipe we present is 
not tied to a particular realization, and could be implemented 
in photonic systems, microwave superconducting circuits, and 
optomechanical systems. 

Finally, while our focus here has been on bilinear interac¬ 
tions between two principle cavity modes, a similar approach 
of balancing coherent and dissipative interactions could be 
used to make nonlinear interactions directional, and could be 
used in more complex cavity lattice structures. Understanding 
how this form of reciprocity breaking leads to useful function¬ 
alities and possibly new photonic states in such systems will 
be the subject of future work. 
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Appendix A: General approach to directionality 


We show here that our basic recipe of balancing coherent 
and dissipative interactions can make any factorizable interac¬ 
tion between two quantum systems directional. Consider two 
bosonic systems 1 and 2, and consider a general interaction 
Hamiltonian of the form: 

■Hcoh (oi02 + ol6[^ . (Al) 

Here, dj is a system j operator (j = 1, 2), implying that 


[At) At)l 







01,0^ 


0{,02 


The operators dj are otherwise arbitrary; in the case where our 
systems are cavity modes, they could be nonlinear combina¬ 
tions of creation/destruction operators, and/or non-Hermitian. 

To make the above general interaction directional, we need 
to introduce its dissipative counterpart. We do this by cou¬ 
pling to a suitably engineered reservoir which couples to both 
subsystems, and which gives rise to a Lindblad master equa¬ 
tion of the form 
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We have a single dissipator with jump operator z = di P 
e^'^ 02 , with a dissipative rate T. The phase ^ appearing in i 
will be kept general for the moment; we see that by tuning it 
correctly, we obtain the desired directionality. 

It is now straightforward to calculate the equations of mo¬ 
tion for some arbitrary operators (n G 1,2) for each sub¬ 
system; as usual, all system-1 operators commute with all 
system-2 operators. Using the above Lindblad master equa¬ 
tion we obtain 
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with (p = cp — The first two terms on the RHS of each 
equation describe the effects of interactions between the two 
systems (both dissipative and coherent), while the respective 
third term describes a local, generalized damping induced by 
the engineered reservoir. To obtain directionality, we take ip = 
Focusing only on the terms coupling the two cavities, the 
EOM take the form: 
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where again, we have dropped the local generalized damping 
terms. This form makes it obvious that for T = A we decouple 
system 2 from system 1: system-2 observables are not influ¬ 
enced at all by system 1, whereas system 1 is influenced by 
system 2. We thus have made the original, general interaction 
described by ?icoh directional, by picking a suitable form for 
the dissipative jump operator and for the corresponding dis¬ 
sipative rate. This simply corresponds to our general recipe 
of balancing a given coherent interaction against its dissipa¬ 
tive counterpart. Note that if we had instead made the choice 
if = — 1 ^ 12 , the direction of the flnal nonreciprocal interac¬ 
tions would be flipped, with cavity 1 now influencing cavity 
2 . 

As an example we consider the directional DPA discussed 
in Sec. |IIC| There we would have the correspondence d\ = 
X 2 and 02 = A • Taking these together with the equations for 
the expectation values in Eq. ( |A5| ), we recover that the QND 
observables X 2 and A are not affected by the interaction, but 
the expectation values for the remaining quadratures, i.e., set¬ 
ting Ai = A and A 2 = Xi, become 

|(p,) = -[A + rl(A). |{.t,) = + |A-rl(x,), 

(A6) 

which coincides with our former result, cf. Eq. 

Finally, the above construction also applies directly to 
fermionic systems, if one takes the operators oj, Aj to be even 
in creation/destruction operators (implying that cavity-1 and 
cavity-2 operators commute with one another). 


Appendix B: Waveguide as an engineered reservoir 

We show here how a simple ID transmission line or waveg¬ 
uide can be used as the engineered reservoir needed in the 
dissipative isolator scheme of Sec. |II A| (i.e., it generates the 
dissipator in the master equation of Eq. ^ with jump opera¬ 
tor z = di + (i 2 )- We take a standard approach for the case 
where the waveguide dispersion is linear over all frequencies 
of interest (see, e.g.. Appendix A of Ref. (471 ). Working in an 
interaction picture at the frequency cCcav of the two principle 
cavities, the Hamiltonian of the waveguide takes the form 

"Hw = fiVG / dx - ko)cR + - ko)cL^ , 

(Bl) 

where ko = cCcav/^G, and we have omitted the explicit x de¬ 
pendence of the waveguide flelds cr{x), cl{x). 

Using the total Hamiltonian H = Hw + ^sb (where T^sb 
is given in Eq. (p^), we And the equation of motion 

(^dt^dx] cr{xA) = ikoCR{xA) 

V^G J 

_^ p 

-\-i - 6(x — Xj)dj, (B2) 

i=i,2 ^ 

with an analogous equation for cl{xA)- The delta-function 
source term leads to a discontinuity in each waveguide held. 


e.g., 


r ^ 

CR{xj +r],t) = CR{xj - 7],t) +i- —(B3) 
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Introducing input and output R flelds associated with cavity j 
in the natural manner, this takes the form of a standard input- 
output relation: 

r ^ 

^R,out{^j i— dj(t). (B4) 

ZVq 

A similar equation (and deflnition of input and output flelds) 
holds for the L held. 

Using the fact that flelds propagate freely between the cav¬ 
ities and I = X 2 — xi > 0, we have 

CR,in [X2 , w] = e®* [‘‘'1 ’■ Cfl,out [ 2 ^ 1 ,^], (B5) 

CL,in[a:i,w] = e®*[‘‘'l*CL,out[a:2,w], (B6) 


where we have Fourier transformed in time and deflned 
k[!jo] = ko ^ uj/vq- We can Anally substitute these results 
into the Heisenberg equations of motion for the cavity opera¬ 
tors dj . Fourier transforming, they take the form 
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where the noise operators are 

^iM = y/v^ (cfl,in[a:i,w] , (B8) 

6M = w] + CL,in[x2, A ' 


Note that we have only retained terms associated with the cou¬ 
pling to the waveguide, as our goal here is to see the form of 
the dynamics it induces for the cavities. 

Consider the case where for all frequencies of uj of interest, 
uj 1/r, where r = 1/vq the propagation time between 
the cavitie s. In this case, we can omit the effects of non-zero 
r in Eqs. (B7), and replace k[uj] by ko. Note that in general, 
we are only concerned with the cavity dynamics on frequen¬ 
cies that are most comparable to r, hence this approximation 
requires rt 1. Once we make this Markovian approxi¬ 
mation, Eqs. ( |B7|) b ecome local in time, and have the same 
structure as Eq. (|10|) and Eq. ( [ST] ) in the main text. Comparing 
the form of the equations, we see that the cavity induces both 
a coherent hopping interaction between the cavities (ampli¬ 
tude Jind = r sin/co//2) and an induced dissipative hopping 
interaction (strength Find = T cos kol). 

The induced coherent interaction here always corresponds 
to a real hopping Jind. As such, we cannot use the waveguide 
to provide both the interactions needed for our isolator scheme 
(i.e., one cannot satisfy the directionality condition of Eq. 
using the waveguide alone). Instead, we can use the waveg¬ 
uide solely to provide the dissipative interaction needed for 
the scheme. We thus require the distance between the cavities 
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to be chosen such that there is no induced coherent hopping 
interaction, i.e., 


kol = riTT, n G Z. 


(BIO) 


If the inte ger n in Eq. ( |B10| ) is even, then the noise opera¬ 
tors in Eqs. (B7) are identical: = ^ 2 - For this choice (and 

in the Markovian limit), the dissipative interactions induced 
by the waveguide are completely described by the dissipator 


jC [ 2 ;] in Eq. ^ with the choice z = di ^ d 2 (as can be shown 
using standard techniques (351). The dissipative interactions 
generated by the waveguide in this limit are thus also equiva¬ 
lent to those generated by the auxiliary cavity implementation 
(cf. Eq. ( p^ ) discussed in Sec. |III AT] Eor n odd in Eq. ( |B10| ), 
one generates the dissipator C[z\ with z = di—d 2 \ upon mak¬ 
ing a gauge change d 2 —^ 2 , this is of course equivalent to 
having z = di d 2 . 
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